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Abstract — In this paper we investigate the problem of allocating 
spectrum among radio nodes under SINR requirements. This 
problem is of special interest in dynamic spectrum access 
networks where topology and spectral resources differ with time 
and location. The problem is to determine the number of radio 
nodes that can transmit simultaneously while still achieving their 
SINR requirements and then decide which channels these nodes 
should transmit on. Previous work have shown how this can be 
done for a large spectrum pool where nodes allocate multiple 
channels from that pool which renders a linear programming 
approach feasible when the pool is large enough. In this paper 
we extend their work by considering arbitrary individual pool 
sizes and allow nodes to only transmit on one channel. Due to 
the accumulative nature of interference this problem is a non- 
convex integer problem which is NP-hard. However, we introduce 
a constraint transformation that transforms the problem to a 
binary quadratic constraint problem. Although this problem 
is still NP-hard, well known heuristic algorithms for solving 
this problem are known in the literature. We implement a 
heuristic algorithm based on Lagrange relaxation which bounds 
the solution value of the heuristic to the optimal value of the 
constraint transformed problem. Simulation results show that 
this approach provides solutions within an average gap of 10% 
of solutions obtained by a genetic algorithm for the original non- 
convex integer problem. 

I. Introduction 

An important aspect of cognitive radio and dynamic spec- 
trum access is how to distribute resources among users. 
Resource allocation in such systems differ from conventional 
wireless systems because the resources available changes 
with time and location, yielding an optimal allocation 
hard. In dynamic spectrum access scenarios, management 
of spectrum is usually divided between two approaches [l]: 
(i) property rights/licensed/exclusive use and (ii) open ac- 
cess/unlicensed/spectrum commons. In the latter all devices 
are treated equal and sharing of the spectrum is done through 
protocols and etiquettes. Although these protocols and eti- 
quettes may be designed to achieve some overall system 
performance, quality of service can not be guaranteed in this 
scenario as everything will be best effort. 

In the licensed scenario, devices can buy and trade access to 
spectrum. In such a scenario there exists central entities which 
grant access to spectrum on a legal basis and central control. 
This approach has led to numerous publications on spectrum 
trading, including demand |2|, pricing |3| and auctions of 
spectrum Q. In such a scenario, quality of service can be 
guaranteed, and in fact must be guaranteed in order for the 



market to function. A natural problem for spectrum holders is 
then how to maximize revenue. 

In this paper we assume secondary users can buy access to 
spectrum for a given time period. When requesting spectrum, 
the secondary users issue a statement about the necessary 
quality-of-service for their operation. In this paper we limit 
our focus to the case where such a QoS metric consists of 
a signal-to-interference plus noise ratio (SINR). Each user is 
willing to pay a certain amount for its service. The goal of the 
spectrum holder is then to maximize the number of users that 
can be granted access to spectrum, while still guaranteeing the 
QoS for these users. This becomes an admission control and 
spectrum allocation problem, which is difficult to solve. 

The difficulty of such a problem lies in the nature of the 
wireless channel: interference. The accumulative nature of 
interference causes the coupling between decision variables 
to be non-linear and usually the optimization problem is non- 
convex (and NP-hard). A central approach is usually assumed 
to solve such problems close to optimality, and techniques such 
as genetic algorithms or geometric programming can be used 
||5J. Another approach is to solve the problem by a distributive 
approach by investigating the Karush-Kuhn-Tucker (KKT) 
necessary conditions for optimality |6|. However, finding the 
KKT multipliers is complex, and in addition the objective 
function must be continuously differentiable and the problem 
has to be convex for the KKT conditions to also be sufficient 
conditions for optimality Q. 

Based on a recent work 0, we use a constraint transfor- 
mation to decouple the joint admission control and spectrum 
allocation problem into two disjoint problems. We show that 
the admission control problem must be solved by a central 
entity in order to guarantee any QoS to the secondary users, 
while the spectrum allocation problem can be solved in a 
distributed manner close to, and in some cases to, optimality. 
This has the advantage of decreasing complexity of the central 
controller and the option of delegating part of the processing 
to the users while still having a central control. 

The rest of this paper is organized as follows: Section |ll] 
provides an overview of related work. Section III defines the 



system model and problem statement. Section IV presents 
the constraint transformation used in access control. Section 
[V] describes the the heuristic algorithm presented in this 
paper. Section [Vl] describes the channel selection algorithm, 
SectionjVnjprovides performance evaluations and Section VIII 



concludes the paper. 

II. Related Work 

Due to the reuse factor used in early mobile phone systems 
and in many of today's broadcast systems, the spectrum 
allocation problem (SAP) has been extensively studied in the 
literature, see for instance S and references therein. However, 
the SAP problem described for cellular networks with reuse 
factors differ from ours by the fact that these problems 
dealt with fixed topology and fixed resource allocations. E.g. 
each base station was allocated x amount of spectrum, the 
problem was finding this x for each base station. Also, for 
cellular networks with fixed topology, solution time is not 
of the essence as the topology does not change. In more 
dynamic networks, such as dynamic spectrum access networks, 
a popular approach to simplifying the problem is to simplify 
the interference model. One frequently used approach is to 
approximate the interference model by a graph with pairwise 
constraints. E.g. node a cannot transmit on the same frequency 
as node b. This is done in for instance fTOl and (TT). However, 
the inefficiency of such graph-based models has also been 
analyzed quite extensively lfT2llfT3l . 

To model interference and successful reception we use the 
physical SINR model |14|, which accounts for accumulated 
interference as opposed to graph-based models. In recent years 
this model has been extensively studied for link scheduling in 
wireless networks, where approximation guarantees on the op- 
timal solution have been shown to exist fT5llfT6llfT7ll . The best 
known solutions yield a constant factor 0(1) approximation 
on optimality. We extend these works by considering multiple 
channels and individual SINR requirements. 

in considers the case of distributing spectrum among users 
with cumulative interference and a universal SINR constraint 
for all users. They assume a common set of channels K 
available to each user where \K\ ^ 1000 and the goal is 
to allocate a fraction of the total number of channels to 
each user according to some objective function. To overcome 
the cumulative interference problem they propose a novel 
linear relaxation of the interference constraints such that the 
problem becomes a linear program. The linear relaxation of the 
interference constraints introduced in fE\ forms the basis of our 
work. Compared to (HI, we extend their work by considering 
individual channel sets and SINR requirements for each user, 
and also assume the individual channels sets can be arbitrarily 
small. Instead of allocating a fraction of a common set of 
channels to each user, we assume a user can obtain a maximum 
of one channel. These assumptions changes the problem from 
a linear program in 1 8 1 to a binary quadratic constraint (BQC) 
problem. 

In 1 18 1 spectrum allocation among femtocells is considered, 
where the goal is to maximize overall system rate while 
guaranteeing a minimum SINR at each femtocell user This 
may seem as a more generalized approach then the one 
considered in this paper However, such a formulation has 
certain drawbacks: (i) with such a formulation each femtocell 
user must transmit, otherwise this problem is equivalent to the 



general sum-rate problem which is known to be unfair (i.e. 
most of the resources are allocated to users with high channel 
gain), (ii) It is known that in many cases some users must stay 
silent for a feasible solution to exist. Thus the requirement 
that each user must transmit on some channel limits the cases 
where a feasible solution exist. However, from a spectrum 
holders perspective, if there is demand then a feasible solution 
exists as one could simply allocate all spectrum to the highest 
bidder 

The work most closely related to ours is ||T9| . In ||T9| 
the goal is to allocate spectrum to as many cognitive users 
as possible while achieving a SINR requirement at each 
user and satisfying some interference hmit at the primary 
users. It is shown that the general problem formulation is 
NP-hard and thus a greedy heuristic is proposed. As the 
problem formulation is still NP-hard in this paper, we also 
rely on a heuristic algorithm to solve the optimization problem 
efficiently. However, we consider a case where primary users 
are only protected through a maximum power constraint at 
each user, such as in the TV white spaces. By this reduction 
we can use a constraint transformation such that we are able 
to utilize heuristics with a much lower complexity bound and 
by comparing sufficient and necessary conditions we are able 
to upper bound the optimal solution of our problem in certain 
cases. 

III. System Model and Problem Formulation 
A. System Model 

A set of N transmitters wants to allocate spectrum to 
support transmission to their receivers. A pair consisting of a 
transmitter and a desired receiver is denoted as a user These 
users can be one-hop links in ad-hoc networks, or they can 
be access points to terminal links. In latter case, if different 
terminals connected to the same access point uses different 
frequencies, a user has to be defined for each terminal to 
support this scenario. However, in general this does not affect 
the spectrum allocation framework, as this only depends on 
how interference between these users are defined. 

To acquire access to spectrum, the users contacts a spectrum 
holder through some means of communication. The accepted 
approach to grant access to secondary spectrum is through 
a database approach as currently done by the FCC [20|. 
We assume the system consists of a set of channels K.. By 
providing location information to a spectrum holder database, 
each user has potential access to a subset ICi of channels along 
with a maximum power constraint P, which guarantees that 
primary users are not affected. Thus we have that ICi C K, 
and Ki — |/Cj|. Each user can select one channel among the 
Ki channels, where each channel has the same bandwidth and 
propagation characteristics. The SINR of user i on channel k 
is given as 



SINRf = 



(1) 



where € {0,1} is 1 if user j transmits on channel k, gj^i 
is the channel gain between node j and node i and cr^ is the 



noise variance. 



Let Si — cji^iPi and Ij i = gj^iPj- The SINR*"' is then given 



as 



SINRf = 



(2) 



We assume each has a SINR requirement /3i. 
Definition 1. Let A define a spectrum allocation, 

A = {a^}, I <i < N,l < k < K, 



(3) 



where ~ 1 indicates that user i selects channel k, = 
otherwise. 

Definition 2. With the SINR requirements of each user, a 
spectrum allocation A is successful if = 1 implies 
SINR'y > A. 

B. Problem Formulation 

As mentioned above each user has a SINR requirement. 
By contacting spectrum holder database, each user is willing 
to pay a fee for access to the spectrum, given that its SINR 
requirement is met. To provide a price function versus time and 
performance is out-of-scope of this paper, we only assume it is 
a non-decreasing function of SINR. The goal of the spectrum 
holder is thus to maximize its revenue by allowing an optimal 
set of users access to the spectrum under their individual SINR 
constraints. 

Since a user can only use one frequency, we introduce a 
variable Xi E [0, 1] where 



The problem can then be formally defined as follows 



s.t.: 



a'liSINR^i - p,) > 



0<x,< l,af G {0,1} 



(4) 

(5) 

(6) 
(7) 

(8) 



where r; is the revenue of the spectrum holder given that user 
i is granted access to the spectrum. This problem is a non- 
convex 0,1 problem, as (|6]l is a non-convex constraint and the 
values of is either or 1 . The difficulty of the problem lies 
in the large unknown non-convex solution space due to the 
large number of possible variable configurations. For instance 
if there are 10 users with 5 available channels each, there are 
5^" possible combinations. 

As noted in numerous papers, this problem belong to the 
class of NP-hard problems. Therefore the time required to 
solve the problem to optimality increases exponentially with 
the number of constraints. The approached used in this paper 
consists of a decoupling of the problem into an access control 
problem and channel allocation problem. The reason for this 



approach is twofold: first, through a constraint transformation 
a sufficient condition for a feasible solution can be stated 
which significantly reduces the solution space and for which 
efficient heuristics exist. Secondly, given a set of users which 
is known to contain a feasible solution, a distributed approach 
to channel selection can be done which in general perform 
close to optimal and in certain cases is optimal. 

Thus the channel selection part of the problem can be done 
both centrally and distributively. This has the advantage of 
being able to offload some of the processing requirements at 
the central entity to the users in the field. On the other hand, 
access control cannot be done in a distributed manner while 
still guaranteeing QoS at the users as "global" information 
about the users is required. In Section |VII| we show how 
limiting knowledge to consist of only certain neighboring users 
affects performance as would be the case in a distributed 
implementation. 

IV. Access Control 
A. Constraint Transformation with Equal Channel Sets 

To overcome the large solution space in the original problem 
(|5]l-(|8]l, we transform constraint (|6| to a binary quadratic 
constraint, which although is not linear, reduces the possi- 
ble combinations of variables to 2^ for which an efficient 
heuristic algorithm can find good solutions. The constraint 
transformation was introduced in IS), but we have modified 
it to our problem. Instead of transforming (|6| from a non- 
convex constraint to a linear constraint, which is possible in 
lis] due to the fractional channel set allocation which renders a 
linear approximation possible, we transform the constraint to 
a binary quadratic constraint due to our single-tone approach. 

The constraint given in (j6]) can be approximated by 



Pi 



Xi{l 

where /™ 

problem then becomes 

max riXi 

i 

S.t. 



0-2 and /+. = min(/™'''^,/j. 



(9) 



The 



(10) 



(1 + Y.^oJ^) < G {0, 1}* - 1, ...,N (11) 



If Xi equals 1, this means user i can achieve its SINR 
by selecting some channel. It does not, however, find the 
channel that should be selected. Let z' be the solution to the 
optimization problem given in ([TOl)-([TT|), i.e. the number of 
users that can obtain their SINR requirement. We then have 
the following proposition: 

Proposition 1. z' is a feasible solution to the original opti- 
mization problem ([5])-([S]). 

Proof See Appendix [A] ■ 



B. Constraint Transformation With Unequal Channel Sets 



V. A Heuristic using Lagrange Relaxation 



Although the solution to ( 10i-( 11 1 finds a feasible solution 



to the original problem, we conjecture that the constraint given 
in (j9]l is too restrictive in some cases. This is especially the 
case when the different users have different subsets of channels 
available to them, which is the case in many secondary 
scenarios as spectral resources depends can depend on location 
and device specifications. In the proof of Proposition [T] a 
feasible solution is shown to exist by bounding the number 
of blocked channels, where all users are assumed to be able 
to block all channels at any given user if they are allowed 
to transmit. However, if two users don't share any common 
channels they cannot block each other. Also, if two users 
only share a subset of channels it is less likely that they 
will contribute the blocking of channels. Thus we propose a 
new constraint for this scenario which takes into account the 
probability of two users blocking each other based on their 
respective channel sets which is yields a feasible solution in 
the asymptotic regime almost surely. The new constraint is 
given by 



/+|/c,n/c 



(12) 



and thus the new optimization problem becomes 



To simplify notation we set — jmaxj^-. — (or — 



if solving (10 1-( 11 1). A Lagrange relaxation of our problem is 

N N N 

maxj^riXi + '^\{Ki - Xi{l + '^a^jXj))] (16) 
s.t.: 



X g {0, 1}^ and A > 0. 



(17) 



Let x' be a solution to (13l-(15l (or (lOi-(lli) with z(x') as 
its value. Let RL(x') be the value of the Lagrange relaxation 
problem. Clearly RL(x') > ^(x'), since when x' satisfies ( 12 1, 
Ki — Xi{l + O'ijXj) ^ 0- Therefore the optimal value of 
( 16 i-lfTT]) is larger than than z* and it is therefore a relaxation. 

It is easy to see that this maximization problem is equivalent 
to the following 



N N 

max{y^(r, - Ai(l + ^aijXj))xi} 

1=1 j^i 

s.t.: 

xe{0,l}^ 



(18) 



(19) 



since XiKi are constants. It is also easy to see that (18l-(19l 
has the trivial solution 



max > TiX 



s.t. 



) < K,,i^ l,...,N 



X, e{o,i},i = i,...,7V 



(13) 

(14) 
(15) 



Proposition 2. Let z* be the solution to the problem given in 
( |7i| )-( |i5| ). z* is a feasible solution to the original optimization 
problem almost surely when N oo. 

Proof: See Appendix [B] ■ 
We now have two new 0, 1 non-linear problems to compute 
the maximum revenue that a spectrum leaser can achieve 
while satisfying the requirements of a successful alloation. 
Specifically, the problems are modified versions of the multi- 
dimensional 0, 1 knapsack problem (MKP), with the difference 
that instead of linear constraints we have binary quadratic 
constraints. As MKPs are NP-hard 121], so are our problems. 
Thus solving the above problems is not easy. By exploiting 
the structure of our particular problems we use a Lagrange 
relaxation inspired by the one introduced in |22| for the 
standard MKP. A heuristic based on Lagrange relaxation has 
desirable properties such as low complexity (running time is 
0{2N^) for our problem) and an upper bound for z* without 
resorting to an LP relaxation of the problem. 



1 if (r, -A,:(l + ^f,^a,,x,)) >0 







otherwise 



(20) 



The difficulty lies in finding values for the Lagrange multipli- 
ers such that the solution x is feasible for the original problem 
and for which the inequality 



N N 

ClijXj ))>0 



(21) 



i=i 



is as close to equality as possible. When equality holds the 
solution X is the optimal solution to the BQC problem [23j|. 

A. An Upper Bound 

Before we start with the process of finding the Lagrange 
multipliers, we provide a result on an upper bound for the 
optimal value of the original MKP. 

Proposition 3. Let x be the solution and A the multiplier 
values of a Lagrange relaxation. Then an upper bound for the 
optimal value of the BQC problem is given by 



N 



N 



z* < z(x) + V X,{K, - x,{l + V a,jXj)) (22) 



i=l 



Proof of this proposition can be found in 11221 Thm 4.1]. 



Algorithm 1 Channel Selection 



1: Consider any x that is a solution to ( 13 1-( 15 1. 

2: For each user i with Xi = 1 select a random channel from 

for each user i with Xi = 1 do 

k ] 



TABLE I 

Set of SINR targets and their revenue value for the two sub 
problems, (i) maximizing number of satisfied users and (ii) 
maximizing revenue 



Compute = apj^,. 
Select k* — arg min^ ijj^ 
end for 

Repeat 3-6 until no more adjustments can be performed. 



B. Finding the Lagrange Multipliers 

Clearly the quality of the solution depends on the Lagrange 



multipliers, as a low value of the summation in ( 22 1 provides a 
small gap to the optimal value. Finding Lagrange multipliers 
for integer programming have been investigated in e.g. [23l 
and |I24|, and specifically for the MKP in |l22l. We adopt the 
algorithm proposed in ll22l . by changing it to our modified 
version of the MKP. 

The algorithm consists of the following steps: 
Step (intialize and normalize) 
Let \i — and xi — 1, i = 1, ...,N. 
Normalize the coefficients: 

Qij —aij/Kj for j = 1, N and i = 1, N 
an =l/Ki for i = 1, ...,N 
b, =1 for i ^ 1, ...,7V 

Compute yi = Ylj'=i for i = 1, N. 

Step 1 (determine the most violated constraint) 

If < 1 for all i, then x is feasible and we can stop. Else 

find the most violated constraint: 

i* — argmaxj/i 

i 

Step 2 (compute the increase of the multipliers Xji , Xji) 
Compute 

j* — argmax ^-^Xo 

Step 3 (increase the Xs) 
Set Xj* = ^j* ^~ di* 

Set Xj^ = 0, y.j* = and yi ^ yi - Uij*. 

If < 1 for all i, then go to step 4; otherwise go to step 1. 

The idea of the algorithm is as follows: find the most 
violated constraint and find the variable that contributes the 
most to violating this constraint as a fraction of the revenue 
it would bring the spectrum holder. Set this variable to zero 
along with its constraint, and reduce all other constraints by 
the contribution this variable had on these constraints. 

VI. Channel Selection 

Through the heuristic algorithm presented in the previous 
two sections we have found the users which are allowed 
to transmit in order to achieve a feasible solution of users 
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that can transmit simultaneously while achieving their SINR 
requirements. However, the solution to this problem does 
not actually find the channels to be used by each user with 
= 1. In a single-channel network where power control is 
used to achieve SINR targets at different users, if a feasible 
solution exists a simple greedy distributed algorithm converges 
|25j. Unfortunately, such a result does not exist for a system 
consisting of multiple channels. 

In the general case this is not a trivial problem. However, by 
assuming that the channel gains between users are reciprocal 
(such that gi j ~ gj.i)^ a simple greedy algorithm is guaranteed 
to converge. Such an algorithm is given in Algorithm [T] 

Proposition 4. Algorithm [7] will converge. 

Proof: See Appendix [C] ■ 
An example of a scenario where such an assumption can 
hold is wireless APs such as Wi-Fi networks. In such a 
network, the APs have about the same coverage area and thus 
the channel gain between AP i and AP j is approximately the 
same as between AP j and AP i. 

VII. Performance Evaluation 

Through this paper we have presented a heuristic approach 
to solving the simultaneous SINR problem. This process 
consisted of two main steps: simplifying the constraints in 
the original problem to binary quadratic constraints (BQCs) 
and using a heuristic algorithm for solving the BQCs. We are 
thus interested in evaluating how the transformation from the 
original non-convex constraints to BQCs affects the optimal 
value as well as how the heuristic performs compared to the 
optimal value of the BQCs. Also, Proposition [T] says that the 
solution obtained from the maximization problem (13l-(15l 
is a feasible solution almost surely for N —i' oo. Thus, for 
finite values of N we are interested in how many of the users 
with Xi = 1 actually achieve their SINR requirements when 
Algorithm [T] terminates. These issues are investigated in this 
section as well as the impact of a distributed implementation 
of the access control algorithm. 

A. Analytical Analysis of Geometric Signal Propagation 

We start by a theoretical analysis of the transformation from 
the problem in Q-dHJ to the one given in (10l-(lll. The 
constraint given in pVis a sufficient condition for a feasible 
solution. If channel gain can be modeled as proportional to 
distance so that the triangular inequalities holcrl and all users 

'let X, y, z be points in a Euclidean space, then d{x, y) < d{x, z)+d{y, z) 
hold. 
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Fig. 1. Solution of tlie original problem (GA) compared to the GA solution 
value of the problem given in (13)-(15) as well as the solution value of the 
BQC problem (lO)-(ll) (GA, heuristic, upper bound and channel allocation) 
for the max sat problem. 
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Fig. 2. Solution of the original problem (GA) compared to the solution 
value of the BQC problem (13)-(15) (GA, heuristic, upper bound and channel 
allocation) for the max sat problem. 



have the same SINR requirement, we can through a necessary 
condition bound the optimal solution as follows. 

Proposition 5. Assume the optimal solution to can 
achieve OPT. Then the optimal solution to (lO)-(ll), OPT', 
can achieve 

OPT 

' (23) 



- - 1 < (9 PT' 
min{2« - 1, 10} 

where a is the pathloss exponent. 

Proof: See Appendix [P] ■ 
Thus, under these conditions we can find a constant factor 
approximation of the optimal solution. However, if received 
signal power don't satisfy the triangular inequalities it has been 
shown that it is NP-hard to approximate the optimal solution 
of a scheduUng with one channel to within iV^^'' for > 
Theorem 6.1]. As adding channels increases the decision 
complexity, this also holds for our problem. 

B. Simulation Setup 

To evaluate the issues mentioned above we simulate an 
environment of constant user density. Specifically, we have 
a density of 1/800 [users/m^]. The rest of the user parameters 



Fig. 3. Actually satisfied users after Algorithm 1 terminates vs solution value 
of the heuristic. 



are set so that the optimal value of the original problem is 
less than the number of users for most instances. We create 
the environment by positioning the users randomly in a square 
according to a uniform distribution in two dimensions. The 
distance from a user to its receiver is given by Gaussian 
distribution with mean value 10 and variance 5. The transmit 
power of each user is set to 1 W and the noise variance is set 
to 10-^ W. 

The SINR requirement of a specific user is chosen at random 
from a set of SINR targets, which is given in Table |l] We 
consider two sub problems, (i) maximizing number of satisfied 
users, in which case the revenue of each user is equal, and 
(ii) maximizing revenue, in which case satisfying larger SINR 
targets leads to a larger revenue. 

To provide a reference value to the optimal value of the 
simplified problem compared to the original optimization 
problem, we solve (|5]l-([8]) using a genetic algorithm (GA) 
from the MATLAB Global Optimization Toolbox. GAs are 
particularly suitable when the structure of the solution space 
of the problem is unknown and little theoretic analysis of 
the problem has been done. To benchmark the quality of the 
heuristic algorithm, we also solve the constraint transformed 
problems through the GA function in MATLAB. The popula- 
tion size of the GA was set to ten times the number of decision 
variables. 

C. Simulation Results 

Fig. [T] shows the solution values for an environment with 5 
available channels where the goal is to maximize the number 
of satisfied users. Since all users have the same spectral 
resources, the constraint transformation in ( [T2] i will not hold 
and as can be seen from the results the solution values obtained 
from ([T3]l-([T5]l are not valid solutions. On the other hand, the 
solutions obtained from the constraint transformation in (9 
are valid solutions, but we can see that the constraint in (9 
is only a sufficient condition for a valid solution as the gap 
between the GA performance of the original problem in (j5]l- 
(|8]l and the performance of ([TOl)-([TT]i is roughly 2 for high 
numbers of users. 

For the results in Fig |2] - 15] we let the spectral resources 
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Fig. 4. Solution of tlie original problem (GA) compared to the solution 
value of the BQC problem (13)-(15) (GA, heuristic, upper bound and channel 
allocation) for the max revenue problem. 
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Fig. 5. Distributed performance as percentage of centralized solution after 
Algorithm 1 terminates. 



at each user vary such that its available channels is drawn 
as a random subset of Ki channels from the total number of 
channels K. In the simulations Ki is a uniformly distributed 
number between 2 and K and K is set to 10. In Fig. |2]the 
solution value of the different algorithms are plotted for the 
max sat problem. For 2 and 3 users, the BQCs optimal value 
can actually be greater than that of the GA because when the 
number of channels at any user is at least 2, the BQC problem 
is trivial for 2 users and will always set a; = 1, as can be seen 
from ( 12 1. For large values of N we see that the solutions 



obtained by BQC problem can surpass the GA solution of the 
original problem. A possible explanation for this is that with 
increased number of users the solution space for the original 
problem is so large that it becomes less likely that the GA will 
find the global optimum. With this in mind, the results indicate 
that the constraint transformation captures the performance of 
the system well. 

In Fig. |3] the difference between the optimal value of the 
heuristic and the number of users actually achieving their 
SINR requirements when Algorithm [T] terminates is given, 
normalized against the former. This further validates the 
constraint approximation done in this paper, as for number of 
users greater than 5 there is a less than 2% difference between 
the optimal value of the heuristic and the number of users 
that actually achieve their SINR. Note that in the simulations 
here we have not assumed reciprocal channel gains (as was 
assumed in the proof of Proposition |4|l, and some of the gap 
can be attributed to this fact. 

In the previous results the goal has been to maximize the 
number of satisfied users. In Fig. |4] the goal is to maximize 
the revenue, where each users SINR requirement yields a 
corresponding profit as given in Table [I] Again, the constraint 
transformation performs within the same percentage of the GA 
solution of the original problem, between 10-15% for 10-18 
users. 

In Fig. |5] we investigate the effect of taking only a subset 
of interfering users into account at each user. As our sufficient 
condition is based on bounding the size of different sets 
of interfering users, limiting the number of users that are 



contributing to the interference at each user will increase the 
number of users the heuristic algorithm allows to transmit. 
The question if how this affects the actual channel allocation 
based on the users that are allowed to transmit from constraint 
transformation. In Fig.|5]we have plotted the number of users 
that actually achieve their target SINR after Algorithm [T| 
terminates for different neighbor degrees, as percentage of 
the number of users that achieve their SINR targets after 
Algorithm[T]terminates with global knowledge at the heuristic. 
Based on the user density we have calculated the distance 
within which x neighboring users are expected to be found. 
Thus, for Neighbors = x, on average x neighbors will be 
taken into account when solving the optimization problem. 
As can be seen from Fig. |5] for a low number of users 
limiting the number of neighbors taken into account actually 
increases the number of satisfied users even after Algorithm [T] 
terminates. Thus, we can assume that the sufficient condition 
is a bit restrictive. However, for increasing number of users 
limiting the number of neighbors that are taken into account 
severely degrades performance. This result limits the sufficient 
condition to be used in a distributed manner. 

The result that capture the benefit of the heuristic algorithm 
the most is given in Table |ll] which shows the time taken to 
solve the problem for the different approaches. As the SINR 
requirements of the users change with time, the time it takes 
the algorithm to find a good solution is important. E.g. we see 
that for 12 users it takes the GA algorithm 30 minutes to solve 
the optimization problem. In this time it is likely that some of 
the users transmission parameters have changed, such as the 
SINR requirement and worst case signal strength, rendering 
the solution suboptimal. 

Due to the time required to solve the original problem 
through a GA, simulation results are only provided for a 
maximum of 18 users. However, as our constraint transfor- 
mation allows us to utilize a heuristic algorithm with a time 
complexity 0{2N'^) this approach can be used to provide good 
estimates of the max sat problem and the max revenue problem 
for large-scale networks. 



TABLE II 

Time to solve the problem 



Number of users 


Algorithm 


Time (s) 


2 


GA 
BQC GA 
Heuristic 


1.982 
0.422 
0.0001 


6 


GA 
BQC GA 
MKP Heuristic 


125.7 
1.072 
0.0004 


12 


GA 
BQC GA 
MKP Heuristic 


1832 
3.922 
0.0009 


18 


GA 
BQC GA 
MKP Heuristic 


16813 
7.305 
0.0015 



VIII. Conclusion and Future Work 

In this paper we have presented a heuristic algorithm for 
solving the spectrum allocation problem under SINR require- 
ments. Specifically, we investigated how many users that can 
transmit and achieve their SINR targets simultaneously and 
maximize the profit to a spectrum holder. As this problem is a 
non-convex integer problem, we transformed this problem into 
two binary quadratic constraint problems, one for which the 
solution is guaranteed to be a feasible solution to the original 
optimization problem and one which is a feasible solution 
asymptotically almost surely. To solve the BQC problems, we 
presented a heuristic algorithm based on Lagrange relaxation 
which bounds the solution value of the heuristic to the optimal 
value of the BQC problem. Through simulation results we 
showed that this approach yields solutions on average at a gap 
of 10% from the solutions obtained by a genetic algorithm for 
the original non-linear problem. 

Appendix A 
Proof of Proposition[T] 



We need to show that for each x. 
k such that SINRj^ 

Let uj^ — J^i^i'^'j-^j-i represent the accumulated interfer- 



: 1 there exists a channel 



ence at user i on channel k. From the definition of I^^^, it is 
clear that SINRj' > j3i is equivalent to cjf < If^^. We define 
a channel k as blocked if oj^ > I™^^. To prove the proposition 
we need to prove that \f Xi — 1 (i.e. user i transmits) there 
exists at least one channel that is not blocked. 

We bound the number of blocked channels as follows 



$ ={fc|a;f > IT'"""} 
,fc 



={fc|3j,a*^ = 1 and Ij^i > /, 
{k\j : a'^ 

= $1 U $2 



^}U 



a-- = 1 and J^- < and w.f > 



(24) 



$1 contains the channels blocked by strong interferers, i.e. 
those that can block a band single handedly, whereas $2 
contains those channels that are blocked by accumulative 
interference. We now bound the sizes of $1 and $2- We first 



bound $1. 



Ki 



<f 1 - E E 4 < E 



(25) 



To bound $2 we use the fact that; 

E = E E 4^.^-^* = E E 
^ E ^j'^^j ^^^^ 

Since for any channel A; e $2 > I™^^ we have that 



iT^ < E ^ E 



fcG'i>2 fee* 
Since J™'*'^ is fixed we have that 
1 



$2 < 



E/ ^3,i^3 



(27) 



(28) 



And thus that the number of blocked channels is at most 

|<i>| = |<fi| + |<i>2|< 

7j, 

. V -^x 

= E^^7=i^ (29) 
I+. 

And thus if ^3 1^ Ki ~ 1 there exists at least one 

channel on which user i can achieve its SINR requirement. ■ 

Appendix B 
Proof of Proposition[2] 

Consider the sets $1 and $2 defined in Appendix |A] We 
now bound them again. We bound $1 as 

I'l^E E 4'^ E ^jProb(j,*) 

(30) 

where Prob(j, i) is the probability that j will contribute to 
blocking on any of i's channels. We note that the same can 
be done for $2 so that 



(31) 



fee* 



and 



fc6*s 

SO that 

I $2 



Y ir"" < E ^ E ^..2;,Prob(j,z) (32) 

fee*2 i#iJj,i<7f""' 



< E IJ,^XJProh{j,i) 



(33) 



3^i,Ij.,<IY"" 



The total number of blocked channels can now be given as and 



|$| = Va;,-^^Prob(?, i) - ^^x,- 



J llCjCiJCil 



Clearly as ^ oo 

/+ 



(34) 



y X, ^Prob(j, z) - y X, \J^llU^ ^ (35) 



due to the central limit theorem. Thus it follows that if 
^/maxj^. ' < Ki — 1, there exists at least one channel 
on which user i can achieve its SINR requirement almost 
surely. ■ 



Appendix C 
Proof of Proposition|4] 

We prove that a greedy algorithm such as the one given in 
Algorithm [T] will converge under the assumption of reciprocal 
channel gains (i.e. gij — gji is the pathloss between user i 
and j). 

Consider a utility function w^' for each user i that depends 
on k as follows 



(36) 



Clearly minimizing in Algorithm [T] is equivalent to mini- 
mizing u'l;. Consider a potential function P(A) defined for an 
allocation A as 



N Ki 



(37) 



=1 fc=i 



We prove convergence by showing that the game forms a 
generalized ordinal potential game. Convergence is then guar- 
anteed as all finite generaUzed ordinal potential games have 
a pure strategy equilibrium [27 Corollary 2.2] and all finite 
generalized ordinal potential games have a finite improvement 
path (|22l Lemma 2.3]. 

Now, user i will change its allocation from channel k to k' 
given that 

= E '^9,^P, < E '^9,^P, = ^ (38) 

By multiplying each side by Pi we have 

= E ^a^PlP^ < E "■J9J^PJP^ = (39) 

Assume before user i changed from k to k' that user I transmits 
on k and user m transmits on k' . Since user i no longer 
transmits on k, clearly ujf (after) < utf (before). Define Aui, 
Aui and Aw^ as 

' = E "■J9J^PJP^ - E 49r^P:>P^ > 0' 



Aujn — Um (before i changed) — (after i changed) 

— ^rn^'iiTiPraPi 

Aui corresponds to the users that have gained by user i's 
change from channel k to k' in terms of increased SINR. Aum 
corresponds to the users that have lost by user z's change from 
channel k to k' in terms of decreased SINR. The total change 
over all affected users is thus 

AP(A, A') = Au, + E + E 

= E a^9l^PJP. - E <^9J^PJP^ + E ""Sj^^J^* 
j¥=i jy^i 

-Y,a^9o^PoP^ 



Y,a'^'g,,P,P,]>Q (40) 



Thus, thus if > uf -P(A) > P(A') and this is thus a 
generalized ordinal potential game. ■ 

Appendix D 
Proof of Proposition[5] 

To prove the proposition we first introduce a necessary 
condition for a successful spectrum allocation: 

Lemma 6. Let x be a successful spectrum allocation. Then 
for each Xi the following holds 



i + E^.]i^)^™(2" + i'io)^* 



(41) 



The proof is similar to the proof of Lemma 3 in HI. 



By comparing the necessary condition of (41 1 with the 
sufficient condition in (j9]) we can bound the optimal value. 

Let be denoted as aj^i and min(2" + 1,10) as C. We 
now define two sets: 



Cs = {i\x^ = l,x^(l - 

^71 — {^l-^i — 



;aj,i) < Ki} 



(42) 



- E 



Thus, \Cs\ is the number of users satisfying the sufficient 
constraints and |£„| is the number of users satisfying the 
necessary constraints. Clearly the optimal value (OPT) of (jSj- 
(|8} is bounded as l^^l < OPT < \Cn\. 

Let ej be smallest a^.j in the sum over £„, i.e. = 
minj^Cn '^j.i- Then we have 

\^n\u < 1 + E ^j'^J'i - '-'^i 

CK,, 



Aui = u; (before i changed)— (after i changed) = afguPiPi 



\Cr, 



which maximizes Consider 

1 + 

= 1 + + ^j'^Jti 

>1+ ^ + (l-CnV-Csl - l)ei 

Thus we have 

1 + ^ Xjttj^i + - l)ej 

<Ki + {\jCn\jCs\-l)ei<CKi 
Simmiation over all i yields 

'^E ^ E(^^ + (.\'^n\Cs\ - l)eO 

>Y,Ki + {\jCn\jOs\-l)Y,ei 

i i 



(44) 



(45) 



C>1 + (|£„\£,|-1) 



C 



(46) 



Since Lg is a subset of £„, |£„\£s| = |£„| — |£s|, and we 



get 



|r„| <C(|A| + 1) 



(47) 
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